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Abstract 

The  compressible  nature  of  the  source  terms  in  Lighthill’s  acoustic  analogy  can  be  closed.  For 
weakly  compressible  flows,  in  the  absence  of  thermoacoustic  effects,  the  compressibility  of  the 
source  field  is  known  in  terms  of  solenoidal  modes  of  the  vortical  flow  field.  In  such  flows,  the 
square  of  the  fluctuating  Mach  number  is  small  and  this  fact,  coupled  with  the  singular  nature  of 
the  acoustic  problem,  and  the  fact  that  the  phase  speed  of  the  acoustic  sources  is  the  advective 
speed,  is  used  to  formally  close  the  compressible  portion  of  the  fluctuating  Reynolds  stresses.  The 
closure  resolves,  as  expressed  in  Crow’s  1970  paper,  the  inconsistent  incompressible  approximation 
to  Lighthill’s  source  term.  It  is  shown  that  the  incompressible  approximation  to  Lighthill’s  source 
term,  accurate  to  order  of  the  square  of  the  Mach  number,  predicts  an  acoustic  field  accurate  to 
order  Mach  number. 
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1.  Introduction 

Many  investigations  in  aeroacoustics  begin  with  Lighthill’s  (1952,  1954)  acoustic  analogy:  it  is 
derived,  without  approximation,  from  the  equations  of  motion  and  constitutes  an  exact  solution  for 
the  acoustic  field.  The  source  term  in  Lighthiirs  equation  while  exact,  requires  in  any  application 
of  the  acoustic  analogy,  approximation.  Any  lack  of  agreement  between  LighthilPs  theory  and 
experiment  is  then  due  to  approximations  made  when  the  source  term  is  modeled.  Lighthill  (1952, 
1954)  showed  that  an  estimate  of  the  source  term  neglecting  its  compressible  nature,  and  where  the 
density  was  taken  to  be  the  local  mean  density,  constituted  a  useful  approximation.  Also  underlying 
Lighthill’s  analogy  is  the  hypothesis  that  the  flow  field,  which  is  the  source  of  the  acoustic  field, 
is  distinct  from  an  ambient  constant  property  stationary  medium  enclosing  the  flow  field;  in  the 
ambient  medium  Lighthill’s  source  terms  are  assumed  to  be  zero.  Several  issues,  limitations,  and 
ambiguities  regarding  the  utility  of  such  assumptions  have  been  investigated  by  Crow  (1970). 

In  this  article  the  compressible  nature  of  the  source  terms  in  LighthilPs  acoustic  analogy  are  closed. 
That  this  might  be  possible,  for  compact  sources,  is  suggested  by  the  singular  nature  of  the  aeroa- 
coustic  problem:  in  the  near  field  source  region,  the  wave  operator,  rescaled  on  the  length  and 
time  scales  of  the  source  flow,  becomes  elliptic,  Crighton  (1975),  Kambe  (1986).  The  compressible 
aspects  of  the  acoustic  source  are  then  given  by  a  series  of  Poisson  equations  whose  source  terms 
are  given  by  the  solenoidal  of  the  inner  source  field.  A  form  of  Lighthill’s  analogy  accounting  for 
the  compressibility  of  the  vortical  source  field  is  obtained  and  the  incompressible  approximation  to 
the  source  terms  in  Lighthfll’s  acoustic  analogy  can  be  understood  and,  if  necessary,  avoided. 

Problem  statement 

LighthilPs  acoustic  analogy  is  written,  Lighthill  (1952),  as 

PM  ~  coo Pijj  =  Tijtij  =  [  pUiUj  +  (p  —  coo p)  8 ij  \iij  •  (^) 

It  is  convenient,  for  the  present  problem,  to  work  in  terms  of  pressure:  LighthilPs  analogy  can  then 
be  written  as 

Oqg  Pm  P >jj  =  (pW'i'U'j') iij  4“  coo  ( P  ^oo p'iitt  (2) 

The  quantity  (p  —  <?00p),tt  is  responsible  for  thermoacoustic  effects,  Crighton  et  al.  (1992);  it 
will  be  neglected  for  the  class  of  flows  that  are  the  subject  of  this  article.  The  present  analysis  is 
understood  to  be  limited  to  the  class  of  flows  relevant  to  subsonic  aeroacoustics:  high  Reynolds 
number  (turbulent),  weakly  compressible,  compact  flows  with  a  relatively  constant  mean  density. 
In  such  a  flow  it  can  be  shown,  using  the  perturbation  series  derived  below,  that  the  quantity 
(p  -  clap)  makes  contributions  that  are  of  0(M?)  with  respect  to  the  highest  order  terms 
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included  in  the  present  analysis.  Here  Mt  =  uc/c  is  the  fluctuating  Mach  number  for  which  uc 
is  a  characteristic  fluctuating  velocity;  its  square  is  a  measure  of  the  strength  of  the  effects  of 

compressibility. 

In  many  analyses  using  Lighthill’s  analogy  an  incompressible  approximation  of  the  source  terms  is 
used:  pUiUj  ~  PooViVj  where  Vi  represents  the  solenoidal  velocity  field.  Lighthfll’s  acoustic  analogy 
is  then  written  as 

coc ?P>tt  —  Pijj  =  Poo{vjVi),ij  (3) 

and  the  difference  between  the  compressible  and  incompressible  term,  [pUiUj  -  PooV^Vj],^,  is  ne¬ 
glected.  This  article  addresses  the  compressibility  correction  term, 

Tfjrij  =  [pUiUj  -  PooViVj],ij  ,  (4) 

which  can  be  viewed  as  encompassing,  in  the  absence  of  thermoacoustic  effects,  the  compressible 
nature  of  Lighthill’s  source  term.  The  present  focus  will  be  on  understanding  the  1)  scalings 
of  Tfi,  2)  identifying  its  acoustic  nature,  and  3)  its  contribution  to  the  radiated  acoustic  field. 
The  main  result  of  this  article,  however,  is  a  closure  for  the  compressible  correction  term,  as  a 
result  the  full  Lighthill  source  term  ( pmuj )>ij  is  then  known.  As  a  consequence,  the  incompressible 
approximation  to  the  source  term  in  Lighthill’s  acoustic  analogy  can  be  analytically  and  numerically 
better  assessed.  It  is  seen  that  in  many  cases  of  practical  interest  the  incompressible  approximation 
is  sufficient.  When  it  is  not  sufficient  the  contributions  due  to  compressibility  can  now  be  obtained. 


Outline 

In  the  next  section,  §2,  some  background  material  and  a  clarification  of  diverse  issues  are  provided. 
In  §3  single-time  low  Mach  number  expansions  of  the  governing  equations  are  conducted.  The 
consequences  and  shortcomings  of  this  procedure  and  its  reliance  on  an  unjustifiable  imposition 
of  a  Helmholtz  decomposition  are  indicated.  In  §4  a  two-time  low  Mach  number  expansions  is 
used;  this  procedure  distinguishes  acoustic  from  compressible  modes  by  phase  speed  -  as  opposed 
to  the  Helmholtz  decomposition  used  by  Crow  (1970)  and  in  §3.  The  analysis  highlights  important 
Vinpmflt.ii?  features  regarding  the  coupling  between  the  solenoidal  and  the  compressible  strain  field 
which  allows  a  simple  and  straightforward  closure  for  the  compressible  portions  of  the  source  term. 
In  §5  the  nature  and  scaling  of  the  compressible  source  terms  are  investigated. 

2.  Preliminaries 

Some  background  information  useful  for  motivating  the  analysis  and  understanding  the  procedure  is 
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now  given.  As  the  singular  nature  of  the  problem  is  in  fact  what  allows  a  closure  for  the  compressible 
correction,  it  is  described  first. 

The  singular  nature  of  the  aeroacoustical  problem 

In  the  present  closure  for  the  effects  of  compressibility  on  the  source  terms  in  Lighhill’s  acoustic 

analogy  the  singular  nature  of  the  aeroacoustical  problem  is  used  to  advantage.  The  singular  nature 
of  the  acoustics  problem  manifests  itself  in  the  nonuniformity,  Lauvstad  (1970),  of  an  expansion 
procedure;  important  physical  processes  take  place  on  different  length  scales,  Crighton  (1975), 
Kambe  (1986),  Kambe,  Minota  and  Takaoka  (1993).  The  ratio  of  the  two  length  scales  goes  to 
zero  as  the  small  parameter,  Mu  goes  to  zero.  When  the  problem  is  scaled  using  the  length  and 
time  scales  of  the  vortical  flow  field,  [t.J./uc],  the  Lighthill  equation,  cj?p)tt  ~  P,jj  =  {(MiUj),ij 
devolves  to  a  Poisson  equation  for  the  pressure  and  one  obtains  a  set  of  equations  that  are  used  to 
describe  an  incompressible  flow:  poo(vi,t  +  VfcU*,*  )  +  P*  =  0,  Vi,i  =  0  and 

V  p  =  Poo{Vi^j\ij  •  (^) 

The  problem  is  elliptic  in  nature.  On  the  outer  or  acoustic  scales,  [f./Mt,f./uc],  one  obtains  the 
sourceless  wave  equation 

cjp,tt  ~  P,jj  =  0.  (6) 

On  the  outer  scales  of  the  problem  the  hyperbolic  nature  of  the  problem  manifests  itself.  The  nature 
of  the  two  problems  is  clearly  different;  Crighton  (1975)  provides  a  pithy  view.  On  the  inner  scales 
the  Laplacian  is  more  singular  than  the  double  time  derivative.  This  fact  is  exploited  to  obtain 
closure  for  the  compressible  portions  of  the  source  terms  will  then  be  related  to  a  series 

of  Poisson  (elliptic)  equations  whose  source  term  is  given  by  the  solenoidal  field  and  is  therefore 
known.  The  essential  point  is  that  the  compressible  nature  of  the  source  term  is  known  without 
solving  the  acoustic  radiation  problem.  In  the  way  that  the  leading  order  contribution  to  the  sound 
field  is  obtained  from  the  solenoidal  modes,  pooViVj,  so  also  is  the  compressibility  correction,  Tfi. 

Observations  by  Crow 

The  neglect  of  the  compressible  portions  of  the  source  term,  T'j ,?J ,  is  sometimes  justified  on  the 
grounds  that  7%  is  of  £>(Mt2)  with  respect  to  (viVj).  This  implicitly  (and  unjustifiably)  assumes 
that  the  compressible  correction,  is  quadrupole  as  is  (i^Uj),#.  Crow  (1970)  has  examined 

this  idea. 

Crow  (1970)  investigated  the  structure  of  the  Lighthill  equation  as  a  perturbation  series  in  the 
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fluctuating  Mach  number. 

M*  p,tt  ~  Vijj  —  'ifjtij  +  ^ijnj  =  Poo(vivj)iij  +  Tij,ij+...  .  (7) 

The  source  term  2-,#  represents  contributions  to  the  source  term  from  interactions  between  incom¬ 
pressible  and,  in  Crow’s  (1970)  nomenclature,  acoustic  modes  of  the  flow.  Crow  (1970)  makes  the 
point  that  to  treat  the  sound  problem  by  retaining  0(M f)  in  the  wave  operator  while  neglecting 
in  Tij  a  compressible  term  of  0(Mt2),  2$,#  is  inconsistent  2.  Despite  his  demonstration  of  this 
inconsistency,  he  concludes  that  the  incompressible  approximation  to  the  fluctuating  quadrupole 
is  adequate  for  compact  flows  with  small  fluctuating  Mach  number.  This  appears  to  be  due  to  the 
nonuniformity  of  the  problem;  p,tt  is  not  important  in  regions  of  the  flow  where  is  important, 
and,  in  regions  of  the  acoustic  field  where  pM  is  important,  T^ij  is  not  important  -  because  of  this 
behavior  the  effects  of  the  truncation  relative  to  the  the  term  p,tt  is  mitigated.  Crow  (1970)  makes 
an  additional  point:  a  simple  source  in  will  contribute  terms  of  the  same  order  as  those  in  the 
incompressible  quadrupole  (i'jVj)ij .  The  major  issue  then  appears  to  be  whether  the  contribution 
to  the  acoustic  field  due  to  Mt2  2$,#  can  be  neglected  relative  to  this  will  not  true  if  2^ 

contains  a  simple  source.  This  important  question  motivates,  in  part,  this  analysis;  its  resolution 
requires  a  delineation  of  the  compressible  source  terms  which  is  only  possible  once  a  closure  for 
these  terms  is  obtained. 

The  governing  equations 

The  following  equations  are  used  to  describe  the  flow: 


P*  )t  +  ukP*  )k  —  ~P*up)P 

(8) 

P  ui)t  P  uk^i  P  n  ® 

(9) 

pVp  00  =  (1 0*}pooV . 

(10) 

An  asterisk  denotes  the  fact  that  the  superscripted  variable  is  dimensional.  The  last  equation  comes 
from  the  assumption  of  isentropy:  ^  s  =  0  =>  pjy-t  h  =  p  from  which  the  equation  follows  given 
the  ideal  gas  law.  The  enthalpy,  h,  is  given  by  h( 7-I)  =  c2  where  <?  =  7 p*/p*  and  =  -§i+u*k^- 
The  momentum  and  continuity  equations  can  be  combined  to  give  p*  ,tt  ~  P*  >j. j  =  {p’<u)U, 
which  becomes  a  wave  equation  for  p  or  p  if  the  gas  law  is  used  to  eliminate  one  in  favor  of  the 
other.  For  clarity  of  exposition,  the  viscous  terms  are  not  carried;  they  can  be  shown,  for  a  high 
Reynolds  number  flows,  to  be  of  higher  order.  Perturbing  about  a  reference  state,  (poo, Poo),  the 
nondimensional  forms  of  the  pressure  and  density  are  taken  as  p*  =  Poo(l  +p)  and  p*  =  Poo(l  +p)- 
2A  similar  inconsistency  in  Ribner’s  (1962)  dilatational  acoustic  theory  has  been  indicated  in  Ristorcelli  (1997). 
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The  equations,  where  [p-,p\  are  now  dimensionless,  become 


Pit  +  w*plfc  — 

1 

)— 1 

+ 

(ii) 

(1  +  pK,t  +utu*,k]  +  ^P,i  = 

0 

(12) 

P-1P  = 

I7(7  _  l)p2 

(13) 

c2 

uoc 

Pitt  ~  P)jj  — 

7 

[(1  +  P)ui  uj]iij  • 

(14) 

The  velocities  are  still  dimensional.  In  these  equations  =  jPoo/poo •  These  equations  will  be 
used  as  the  starting  point.  The  following  definitions  for  the  strain  and  the  rotation  will  also  be 
needed:  Sij  =  ±[vuj  +vj)i  ]  and  m  =  \[vhj  -vjfi].  The  superscript  “c”  as  in  4  and  4  will  denote 
the  strain  and  the  rotation  associated  with  the  compressible  portion  of  the  velocity  field, 

3.  A  single  time  perturbation  analysis 

A  small  fluctuating  Mach  number  expansion  of  the  compressible  Navier- Stokes  equations  is  now  con¬ 
ducted.  The  analysis  is  understood  to  be  addressing  high  Reynolds  number,  weakly  compressible, 
turbulent  flows  -  the  class  of  flows  relevant  to  subsonic  aeroacoustics  for  which  isentropy  implies 
dp  =  <?dp.  Similarites  to  the  Crow  (1970)  analysis  are  discussed.  Shortcomings  are  highlighted. 

A  single  time  perturbation  analysis 

Regular  expansions  of  the  form  p  =  e2  [  p\  +  e2P2  +  •■•]»  P  =  ^  [  Pi  +  e2P2  +  •■■]>  ui/^c  =  Vi  +  wi  = 
vi  +  e2  [  wu  +  e2^2i  +  ...],  are  inserted  into  the  governing  equations.  To  the  lowest  order  a  set  of 

equations  that  are  the  same  as  that  used  to  describe  an  incompressible  flow, 


Vi,t  +VpViip  +  Pin 

=  0 

(15) 

Vi,i 

=  0 

(16) 

V2pi 

=  —  (ViVj),ij 

(17) 

Pi 

=  7P1 

(18) 

are  obtained.  The  velocity  has  been  made  nondimensional  with  uc  and  e2 

=  7 Mt2  =  7^c/coo  where 

=  7Poo/poo-  The  independent  variables  are  still  dimensional  —  length  and  time  scales  have  not 
been  chosen.  The  fluctuating  pressure,  p\  is  obtained  from  a  Poisson  equation  whose  source  term 
is  known  in  terms  of  the  solenoidal  modes  of  the  velocity  field;  the  leading  order  density  fluctuation 
is  linearly  related  to  the  leading  order  pressure.  The  next  order  equations  are 

Pl,t  +  VpPup  =  -W\k,k  =  ~  di  (19) 
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Wli,t  +  VpWu,p  +WipVi,p  +  P2,i  =  -  pl(Vi,t+VpVi,p)  (20) 

P2~'YP2  =  2  7(7 -IV?  (21) 

-  P2,jj  =  {wuVj  +  wijVi  +  p\ViVj),ij  -  -Phtt  ■  (22) 

The  compressible  portion  of  the  source  flow,  [wu,pi,p2],  does  not  involve  a  wave  equation.  In  the 
near  field  region  the  compressible  pressure  is  felt  effectively  instantaneously.  Since  7 pi  =  pi,  the 
next  order  continuity  equation  yields  a  diagnostic  relationship  for  the  dilatation,  d\  =  w\ k,k, 

-'fdi  =  pi,t  +VkPi,k  =  -^Pi-  (23) 

where  -Rj  =  4  +  vk-£~.  The  dilatation  is,  to  leading  order,  not  an  acoustic  quantity  -  it  does 

Dt  ot  f/Xfc 

not  propagate  with  the  acoustic  phase  speed.  The  dilatation,  d\  =  wu,i  is  related  to  the  pressure 
that  satisfies  V2pi  =  —  It  can  be  shown,  Crow(1970),  that  the  solenoidal  pressure  falls 

off  from  the  source  region  as  pi  ~  x~3.  Using  the  continuity  equation  indicates  that  uq*,*  =  d\  ~ 
■§1  Pi  ~  x~3  and  wu  ~  x~2  which  differs  from  the  acoustic  scaling,  wu  ~  aT*1. 

LighthilPs  acoustic  analogy  is  reconstituted,  in  dimensional  form,  as 

cjp,tt  -  V2P  =  Poo  [ViVj]  ,ij  +  Poo  [WiVj  +  WjVi  +  pViVj]  ,ij  .  (24) 

The  identifications  Wi  =  nfM^ucwn  and  p  =  yMfpi  have  been  made.  This  is  what  Crow  (1970) 
has  called  a  self-consistent  form  of  Lighthill’s  acoustic  analogy  —  0(M?)  terms  on  both  sides  of 
the  equation  are  retained. 

Imposition  of  a  Helmholtz  decomposition 

In  order  to  close  the  source  terms  expressions  for  pi  and  wu  are  required.  This  is  readily  accom¬ 
plished  as  the  leading  order  problem  [v^pi]  is  considered  known. 

In  the  perturbation  anslysis  above  the  small  parameter  was  used  to  order  various  dynamical  aspects 
of  the  physics.  The  problem,  following  Crow (1970),  is  now  kinematically  partitioned:  the  Helmholtz 
decomposition  allows,  for  any  suitably  compact  vector  field,  a  unique  decomposition  into  the  curl 
of  vector,  u*,  and  the  gradient  of  a  scalar,  ui  =  Vi  +  Wi  —  Vi  +  <f> The  decomposition  implies, 
of  course,  wu  =  and  V2<£i  =  d\  and  the  usual  Biot-Savart  relation  V2Ui  =  Sijk^k^j-  Here 
a =  EijkVj.k  is  the  vorticity.  The  compressible  modes  are  therefore  determined  by  the  inversion  of 
the  following  Poisson  equation  for  the  potential 

— 7V2<£i  =  Pi,  (25) 
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using  the  methods  of  Green’s  functions.  The  homogeneous  solution  is  the  trivial  solution.  The 
velocity  field  is  given  by 

-7VVi  =  (§-t  PiU  =*  wu  =  -±v  2(§-t  PiU-  (26) 

There  are  several  other  simplifications  that  can  be  made  to  this  expression  using  the  fact  that  wiX 
is  a  double  convolution.  This  does  not  suit  present  intentions.  The  fluctuating  density  is,  of  course, 
linearly  related  to  the  fluctuating  pressure  jp\  =  p\.  Thus  all  terms  in  [wuVj  +  w\jVi  +  p\ViVj ]  are 
closed  once 


V2P1  =  -  [viVj],ij  (27) 

is  solved.  The  Lighthill  equation,  to  order  M2,  has  now  been  formally  closed.  This  is  similar 
to  Crow’s  (1970)  analysis;  Crow  (1970),  however,  was  interested  in  different  questions  and  did 
not  pursue  the  question  of  a  closure  for  the  compressible  correction  T^,ij  =  [ puiUj  —  = 

Poo[WiVj  +  WjVi  +  pViVj\,ij . 

There  are  some  theoretical  and  practical  problems  with  this  result. 

•  The  Helmholtz  decomposition  has  partitioned  the  flow  field  such  that  all  the  vorticity  is  in  the 
leading  order  problem,  [v*,pi];  all  the  effects  of  compressibility  are  in  the  irrotational  mode. 
The  Helmholtz  decomposition  precludes  any  representation  of  a  compressible  rotational  mode. 

•  Inspection  of  the  evolution  equation  for  wu  shows  that  v>u  is  not  irrotational. 

•  Contributions  to  the  source  term  wu,j  Vj,i  =  (sfj  -~r^j)vj.i  does  not  account  for  the  interaction 

term  since  the  rotational  compressible  are  not  captured  in  this  representation.  Here 

Sij  =  +wlj’i]  and  rij  =  jK-j  ~wljn] 

•  Any  calculation  of  the  compressible  correction  requires  a  double  convolution.  Any  solution 

for  the  acoustic  field  requires  a  triple  convolution. 

These  shortcomings  can  be  resolved  with  a  multi-time  analysis. 

4.  A  two-time  perturbation  analysis 

The  ability  of  a  two-time  expansion  to  distinguish,  in  the  source  field,  compressible  from  acoustic 
modes  is  essential  to  closing  for  the  compressible  aspects  of  the  source  terms  in  Lighthill’s  acoustic 
analogy.  The  multi-scale  procedure  distinguishes  acoustic  modes,  with  acoustic  phase  speed  from 
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compressible  modes,  with  advective  phase  speed.  This  is  a  distinction  that  is  not  possible  with  the 
methodology  employed  above. 

In  the  flow  field  there  axe  two  velocity  scales,  the  sound  speed,  c,  and  a  characteristic  fluctuating 
velocity,  uc.  In  the  flow  field  one  length  scale  is  recognized:  L  Two  time  scales  can  be  identified: 
a  fast  time  that  scales  with  the  eddy  crossing  time,  £/coc,  and  a  slow  time  that  scales  with  the 
eddy  turnover  time,  t/uc.  The  ratio  of  these  two  time  scales  is  the  turbulent  Mach  number  and 
is  the  small  parameter  in  a  two-timing  procedure.  The  problem  is  recognized  as  being  driven  by 
the  vortical  flow  which  evolves  on  the  advective  (slow)  time  scale:  if  there  were  no  flow  evolving 
on  the  slow  time  there  would  be  no  compressible  or  acoustic  field.  The  problem  can  be  viewed  as 
a  forced  linear  oscillator  with  forcing  coming  from  the  slow  vortical  modes  with  advective  scales 
[ucJ].  While  simple  in  concept  diverse  subtelties  associated  with  the  aeroacoustical  problem  can 
be  more  readily  seen  if  the  two-time  procedure  is  first  applied  to  the  forced  linear  oscillator. 

The  forced  linear  oscillator 

Consider  the  forced  linear  oscillator  driven  with  frequency  u  much  slower  than  the  natural  frequency 
of  the  oscillator,  cjo,  <<  1- 

U  +  (Jq  U  =  UqCOS  Lot,  (28) 

The  exact  solution  is  written  as  a  sum  of  homogeneous  and  particular  solutions: 

LJ2 

u(t;cjQ,a;)  =  uh^loq)  +  up^loq^lo)  =  ^()i  H — 2  0  ^  cos  cot,  (29) 

ClJq  LO 

Note  that  the  homogeneous  solution  has  no  relation  whatsoever  to  the  forcing,  uH(t]u o)  /  /(w). 
This  property  of  the  homogeneous  solution  will  be  used  to  eliminate  its  consideration  as  a  solution 
relevant  to  the  acoustic  source  generated  by  the  flow  field.  If  the  initial  conditions  are  such  that 
the  eigensolutions  are  not  stimulated,  A\  =  A2  =  0,  the  solution  is  the  particular  solution,  u(t)  = 

$j^  =  cos(wt)  +  0(fy. 

The  forced  linear  oscillator  is  now  treated  by  the  method  of  multiple-scales.  The  role  a  small 
parameter  will  be  played  by  e  =  w/w o  «  1-  Time  can  be  rescaled,  t'  =  wq/.  to  produce,  after 
dropping  the  primes, 

u  +  u  =  cos  et.  (30) 

Following  the  usual  multi-scale  procedures  the  original  time  variable,  [t],  is  replaced  by  two 
independent  time  variables,  [to,*i];  to  is  the  fast  time  scale  and  t\  is  the  slow  time  scale.  The 
multi-time-scale  ansatz  for  the  dependent  variable  is  u (x,  t)  =  u(x,  ty,  t\),  where  to  =  t  and  t\  =  et. 
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The  time  derivative  of  u(x,  t)  is  then  written,  using  the  chain  rule,  ^  u(t;  e)  =  u(to,ti',e)  + 

^tW1u(t o,h',e)  =  ^u(t0)tue)  +  e-^u(t0,tue).  Expanding  u(t;e)  =  u0(t0,ti)  +  eui(t0,ti)  + 
e2u2  (to,  ti)  +  ...  the  problem  for  the  forced  oscillator  takes  the  form 

■^pUo  +  ua  =  costi  (31) 

w0u'+u'  =  ~2Bmina-  (32) 

The  leading  order  solution  can  be  written  as  a  sum  of  homogeneous  and  particular  solutions, 

u0(t0,t i)  =  A1(t1)e+itn  +  A2(t1)e  u"  +  cos  h.  (33) 

The  slow  time  coefficients,  Ajjti),  are  determined  by  the  removal  of  secular  terms  in  the  equation 
for  u\.  They  are  constants  for  this  problem.  The  solution,  when  the  eigenmodes  are  not  stimulated, 
A\  =  A?  =  0,  is  the  particular  solution:  v.o(to,  ii)  =  Uop(tx)  =  costi  —  cosut  which  is,  to  0(e2), 
the  exact  solution  given  above.  The  realization  that  the  particular  solution  u0(ti;  oS)  is  the  solution 
relevant  to  the  forcing  and  that  the  homogeneous  solution,  iio/f(to;  wo),  is  not  associated  with  the 
forcing  and  not  relevant  to  the  aeroacoustical  problem  driven  by  the  fluctuating  velocity  field. 

Multi-scaling  the  compressible  source  flow 

A  multi-scaling  of  the  compressible  Navier  Stokes  equations,  modulo  viscosity,  is  now  conducted. 
An  analogy  between  the  compressible  aspects  of  the  source  flow  and  the  forced  linear  oscillator  is 
made.  The  time  scale  of  the  forcing,  w_1,  can  thought  of  as  an  eddy  turnover  time  of  the  turbulent 
source  field,  lfuc.  The  slow  or  advective  time  scale  is  ijuc.  The  time  scale  u>a  1,  can  be  identified 
with  the  sound  crossing  time,  i/c^.  The  fast  or  acoustic  time  scale  is  f,/uc.  The  small  parameter 
is  the  turbulent  Mach  number  which  is  the  ratio  of  these  two  time  scales. 

Compressible  portions  of  the  flow  evolving  on  the  fast  time  scale  will  be  called  acoustic  modes; 
compressible  portions  of  the  flow  evolving  on  a  slow  time  scale  will  be  called  advective  compressible 
modes.  The  multi-time  procedure  is  crucial  to  distinguishing  these  two  aspects  of  the  compressible 
flow.  The  compressible  portions  of  Lighthill’s  source  terms  will  be  seen  to  be  a  function  of  the 

advective  compressible  modes. 

Starting  with  the  compressible  Navier  Stokes  equations  given  in  §2  the  following  ansatz,  as  suggested 
by  earlier  results,  for  the  dependent  variable  is  proposed: 


U* 

—  Uc(Ui  -f-  tfiw'Mi) 

(34) 

p 

=  4>a{Ps  +  4>ipc) 

(35) 

p 

=  4>a{Ps  +  <j>iPc)- 

(36) 
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The  subscript  V  is  understood  to  indicate  solenoidal  in  reference  to  which  has  zero  divergence 
and  from  which  both  ps  and  ps  are  derived.  The  subscript  “C”  will  signify  the  compressible  nature 
of  the  nondimensional  pressure  and  density  fields.  Inserting  the  decompositions  into  the  governing 
equations  produces  on  the  inner  length  and  time  scales  the  problem: 

Viit  +  VkVj,k  +  ps,i  =  0  Vk,k  =  0  (37) 

Psijj  =  ~  (ViVj),ij  IPs  =  Ps-  (38) 

The  leading  order  problem  is  described  by  equations  used  to  model  an  incompressible  flow  with 

characteristic  time,  length  and  velocity  scales  [^J.,uc}.  The  gauge  function  <f>a  =  iM?  is  established 
by  balance. 


The  compressible  portion  of  the  problem,  \pc>PC>wi\>  is 

[1  +  <po{ps  +  4>lPc)]\wi,t  +VkWi,k  +WkVi,k  +4>wWkWi,k  \4>w  + 

=  ~4>o{Ps  + 

<h4>i[Pc,t  +vkPc,k  +4,wWkPc>k  ]  +  [l  +  <Po  (Ps  +  0i  pc)]wk,k  4>w  = 

<I>o4>i(pc  -  ipc)  = 


—<Po<PiPc,i  = 

7 

4>iPc)[vi,t+vm,k]  (39) 

-4>o[Ps>t+VkPs,k]  (40) 

i7(7  _  i)4>l(ps  +  4>iPC)\a\) 


The  dependent  variables  of  this  problem  are  all  nondimensional:  [vi,Wi,ps,pc,  Ps,Pc]  are  in  nnits 
of  [uc,poo,Poo]-  It  is  important  to  note  that  the  independent  variables  are  dimensional :  length  and 
time  scales  have  not  yet  been  chosen.  The  scalings  for  length  and  time  will  determine  <pi  and  <f>w. 
Choosing  the  £/c  and  t  for  the  scalings  the  only  consistent  balance  of  the  momentum  and  continuity 
equations  produces  <j>w  =  M2  and  <j>\  =  Mt- 


The  compressible  independent  and  dependent  variables  are  now  expanded  in  a  series  of  form: 


m 

=  ptQ  ffato,h)  +  /(®,to,ti) 

(42) 

Wi(x,t0,ti) 

=  woi(x,t0,ti)  +  ewu(x,t0,ti)  +  ... 

(43) 

pc{x,tQ,h) 

=  po(x,to,ti)  +  epifato,ti)  +  ... 

(44) 

Pc  fa  to,h) 

=  po{x,to,t\)  +  epi(x,to,ti)  +  ...  . 

(45) 

Here  rk>  =  dul  and  =  dk  +Vkak'  The  e9uations  are  made  nondimensional  with 

the  integral  length  scale,  £,  and  the  fast  time  scale,  f./c.  Inserting  the  expansions  for  the  dependent 

variables  and  the  temporal  differential  operators  produces  to  leading  order: 


D 

Dto 


Wo  i  +  Po  ,i 


D 

jy-  Pa+wok,k 


P  o  -  7Po 


0 
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0. 


(46) 

(47) 

(48) 


The  set  of  equations  are  equivalent  to  those  of  the  forced  linear  oscillator,  y’i— 3/2  =  0  an(l  J/2+^o  3/i  = 
cos  ( u>t ).  The  analogy  can  be  seen  in  Fourier  space  where  the  complications  of  the  spatial  derivative 
on  wot  are  removed.  (The  second  order  equations  are  given  in  the  appendix.)  The  compressible 
velocity  gradients  evolve  according  to: 


D 

Dt0 


POyij  ^pyj  ^Uiyp  i 


(49) 


from  which  the  equations  for  the  divergence  and  the  vorticity  are  readily  obtained, 


D 

Dto 


WQj  ,j 


D 

Dt0 


^0  q 


Po  )jj  Vp  )j  ^0 j  yp  • 

tqijVpyj  WBiyp  • 


The  associated  wave  equations  are  given  in  the  Appendix. 


(50) 

(51) 


The  particular  solution  and  its  consequences 

Recalling  the  forced  oscillator  analog,  the  general  solution  for  the  leading  order  problem  is  written 
as  a  sum  of  homogeneous  and  particular  solutions,  wa  =  w^h {x>  to)  -\-wup{x,  ti),  po  =  Pqh(%,  to)  + 
Pip(x,  ti).  The  homogeneous  solution  is  independent  of  the  forcing  and  is  relevant  to  acoustic  fields 
that  are  generated  by  initial  or  boundary  conditions  and  is  not  of  interest  to  the  acoustic  source 
problem. 

The  portion  of  the  compressible  problem  driven  by  the  slow  solenoidal  modes  is  of  interest;  the 
solution  of  relevance  is  the  particular  solution,  wqi  =  WQip(x,  ti),po  =  Pop(xy  ^l)-  As  was  seen  in  the 
forced  linear  oscillator  example,  the  particular  solution,  to  0(e 2),  is  found  by  setting  ^  [u;o,Po]  = 
0.  Thus: 

WQkyk  =  -  Pa  (52) 

P0OM1)  =  PoOMi)  =  0-  (53) 

The  two-time  expansion  indicates  slow  compressible  modes  are,  at  this  order,  irrotational,  and 
found  by 

w0k,k=  V2^(.T,fi)  =  Ps-  (54) 

Rotational  modes  can  exist  at  this  order;  they  are  associated  with  the  homogeneous  solution  and 
axe  a  function  of  arbitrary  initial  conditions.  The  equation  for  the  divergence  with,  -j^-  w0 j,j=  0, 
indicates  that  the  particular  solution  has  the  following  property: 

Vp,j  WQj  ,p  =  Vpij  <j>ijp  =  0.  (5b) 
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This  is  a  consequence  of  the  two-time  expansion;  it  is,  as  will  be  seen,  a  very  useful  fact.  Kine¬ 
matically  it  is  understood  as  indicating  that  vp  yj  WQj  ,p  =  spjs<jp  =  0;  the  principal  ’ axes  of  the 
compressible  strain  rate  are  orthogonal  to  the  principle  axes  of  the  solenoidal  strain  field.  In  the 
compressible  portions  of  the  acoustic  source  quantities  such  as  (wuVj)^j  can  be  simplified  using 
the  above  fact:  (vjfa  ),y  =  Vj  V2<£,j  +  v^j  4>rij  =  Ps)>j  smce  vp,j  <f>,jp  =  0.  The 

compressible  correction  term,  Tfj  nj  =  Poo  [wmVj  +  wojVi  +  P\vivj\  rij  —  Poo  [<f>n  vj  +  j  vi  +  Ps^ivj]  nj  > 
is  then  written  as 

Tij  rij  =  Poo  [—  2  (Vj  Ps)  )j  +  (Psvivj)  rij  }  (^®) 

and  LighthiH’s  acoustic  analogy  is  reconstituted,  in  dimensional  variables,  as 

cjp,tt  -  V2P  =  Poo{viVj),ij  +  ^  cJ[{psViVj),ij  -  2 (v,~  Ps),j}-  (57) 

As  the  perturbation  analysis  is  finished,  is  now  dimensional.  The  compressible  nature  of  the 
source  terms  can  be  obtained  from  the  solenoidal  field  and  the  equation  is  closed:  ps  is  obtained 
from  its  Poisson  equation;  -V2ps  =  Poo  (viVj)rir  Several  featuras  and  advantages  of  the  multi-scale 
closure  are  summarized: 

•  The  compressible  field  is,  to  leading  order ,  irrotational.  The  multi-scale  result  justifies  the 
assumption  of  a  Helmholtz  decomposition  in  the  previous  section  (and  in  Crow’s  (1970) 
analysis). 

•  The  principal  axes  of  the  compressible  strain  rate  are  orthogonal  to  the  principle  axes  of  the 
solenoidal  strain  field,  since  vp)j  wqj)P  =  0. 

•  The  compressible  correction  term  can  be  written  in  terms  of  the  solenoidal  modes  of  the 

source  flow: 

Tij  rij  =  ~  Coo  [(PsViVj)rij  — 

•  The  solenoidal  field,  [vi,ps\ ,  is  calculated  from  the  leading  order  equations  —  a  set  of  equations 

that  is  describes  an  incompressible  flow. 

•  The  compressible  correction  term  does  not  require  a  double  convolution. 


Caveat 

The  analysis  presented  is  relevant  to  a  general  fluctuating  flow  in  which  there  is  one  characteristic 
fluctuating  velocity  scale  and  whose  characteristic  Mach  number,  the  turbulent  Mach  number,  is 
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small.  For  such  a  flow  the  slow  compressible  field  is,  to  leading  order,  irrotational.  The  irrotation- 
ality  of  the  slow  compressible  field,  for  subsonic  flows,  is  a  useful  approximation  as  long  as  the  mean 

strain  rate,  S,  does  not  approach  the  eddy  crossing  time:  it.  when  St/c  — *  1.  When  SC/ c  ~  1,  the 
slow  compressible  modes  become  rotational  and  SC/ c  enters  the  analysis  as  an  additional  governing 
parameter.  Thus  the  present  analysis  is  relevant  to  applications  for  which  SC/ c  <  1;  this  is  the  case 

for  many  developed  turbulent  flows  for  which  uc  ~  SC. 

5.  Detailing  the  compressible  contribution  to  the  acoustic  source 

As  the  analytical  work  is  now  complete  the  presentation  returns  to  the  dimensional  variables  as 
used  in  §1.  The  main  result  of  this  work  has  been  the  closure  for  the  compressible  correction  to 
the  Lighthill  source  terms, 

1 fj,ij=[pUiUj-  PooViVj^ij  =  -^-[{PsViVj),ij-2(Vj^-ps)>j]. 

7coo  Ul 

Various  properties  of  this  closure  for  the  compressible  source  terms  are  now  investigated.  The  far 

field  solution  to  the  wave  equation  is  presented  and  the  nature  of  the  new  source  terms  is  defined. 
The  asymptotic  scaling  of  each  of  the  terms  is  also  indicated.  It  is  first  noted  that  the  second 

compressible  source  term  can  be  rewritten 

(vjd),j  (59) 

as  the  flux  of  the  leading  order  dilatation  using  the  diagnostic  relation,  d  =  -  ^  ps,  derived  above. 
Note  that  if  the  pressure,  ps,  were  frozen,  it  would  not  contribute  to  the  compressible  source  term. 


The  nature  of  the  compressible  contribution  to  the  source 

The  contribution  of  the  compressible  source  terms  to  the  acoustic  field  is,  however,  determined  by 
integrating  over  the  volume  of  the  source  field.  The  closed  Lighthill  analogy  reads 

cjp,tt  ~  V2P  =  Poo (viVj) ,ij  +  f  cj [(psViVj) ,ij  -  2 (Vj  jytPs),j],  (60) 


and  the  solution,  for  an  unbounded  flow,  is  obtained  by  convolution, 

1  f  dhd 

P(M)  = 

1  1  d2  f,  ,  dyf 

+  47T7  dxidxj  J  PsVlVj  |x  -  x' 

4tt7  ci,  dxj  J  Vj  Dt  Ps  |x  -  x'| 


(61) 
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where  the  square  brackets  indicate  a  quantity  evaluated  at  [x;,  tr]  where  tr  —  t  — 
retarded  time. 


.  is  the 


|x-x'l 

Coo 


Following  the  usual  procedures  for  the  far  field  solution,  due  to  symmetry  in  space 

and  time.  The  fax  field  solution  is, 


p(x,  t)  = 
+ 
+ 


_1 

47r  c2 


1  1  XiXj  d2 


x2  dt2 
2  11  xa  d 


47T7c; 

1 


J  Poo[ViVj]  dx 

J [PsViVj]  dx 


1  1  XiXj  d2 


4w~f  x  x 2  dt2 


(62) 


From  a  mathematical  point  of  view  the  three  contributions  to  the  acoustic  field  —  ordered  according 

to  the  expansion  parameter  —  are  seen  to  scale,  respectively,  as  0(M2)  +  O(Mf)  +  0{Mf ) .  It 
is  seen  that  though  T*  is  0(M 2)  with  respect  to  ViVj ,  in  Light  hill’s  wave  equation,  its  contribution 
to  the  acoustic  field  is  0(Mt )  with  respect  to  the  contribution  from.  ViVj.  From  a  more  applied 
viewpoint  the  contributions  to  the  acoustic  pressure  scale  as 


p(x,t)  ~  0(ui)  +  0{ul)  +  0{ul). 


(63) 


The  first  term  produces  the  well  known  eight  power  velocity  scaling  for  the  sound  intensity  of  a 

quadrupole  source,  Lighthill  (1952).  To  leading  order  (pp)  ~  it®;  to  next  higher  order  crossterms 
will  produce  a  u\  scaling. 


The  second  term,  involving  the  flux  of  the  dilatation,  can  be  understood  as  a  lateral  quadrupole. 
Using  the  momentum  equation  the  flux  of  the  dilatation  can  be  rearranged  as 


VjWt  Ps  =  ^  (VjPs^  +  ’ 


Dt 


(64) 


giving  a  clearer  indication  of  quadrupole  nature.  The  (psps),j  will  make  a  higher  order  contribution, 
at  the  same  order  in  Mt  as  the  last  term.  The  second  term  will,  nonetheless,  be  left  as  a  flux  of 
the  leading  order  dilatation.  The  last  source  term  is  a  quadrupole. 


Asymptotic  decay  of  the  compressible  portion  of  the  source  field 

As  a  consequence  of  the  Biot  Savart  law,  Crow  (1970),  Howe  (1975),  Kambe  (1986),  the  boundedness 
of  the  first  moment  of  the  vorticity  indicates  that  the  fall  off  of  the  solenoidal  velocity  associated  with 
a  compact  vorticity  distribution  is  ~  x  3.  As  a  consequence,  the  incompressible  approximation 
to  the  source  falls  off,  with  distance  from  the  source  region,  as  ViVj  ~  x~6. 

Similarly  it  has  been  shown,  see  Crow  (1970),  that  ps  ~  x  3.  As  a  consequence 

W0j  ,j  ~  Ps~X  3  (65) 
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and  Wi  r*j  M?x  2.  Using  the  above  facts  the  source  terms  in  Lighthills  acoustic  analogy,  ViVj  + 
WiVj  +  WjVi  +  psViVj,  fall  off,  from  the  source  region,  as: 

ViVj  ~ 

WiVj  ~ 

PsViVj  ~ 

Thus  in  the  inner  source  flow  there  are  contributions  due  to  the  compressible  nature  of  the  fluctu¬ 
ations  that  are  both  more  and  less  compact  than  the  incompressible  approximation  to  the  source 
term.  The  assumption  underlying  Lighthill’s  acoustic  analogy,  that  Tg  =  0  outside  the  source 
region,  is  validated  when  the  compressible  nature  of  the  source  term  is  included. 

Note  that  the  compressible  velocity  field  scales  as  a;-2;  this  is  not  the  acoustic  scaling  assumed  in 
Crow  (1970).  Crow’s  (1970)  arguments  were  based  on  the  fact  that  the  compressible  component  of 
the  velocity  scaled  as  aT1  thus  w^j  ~  x~4.  As  has  been  shown  the  compressible  correction  term, 

Tc.  is  due  to  the  advective  nature  of  the  compressible  flow,  not  its  acoustic  nature. 

6.  Summary  and  conclusions 

For  a  high  Reynolds  number  weakly  compressible  flow,  in  the  absence  of  important  irreverssible 
processes,  a  multi-time  scale  procedure  has  been  used  to  distinguish  compressible  modes  with  acous¬ 
tic  phase  speed  from  compressible  modes  with  advective  phase  speed.  The  advective  compressible 
modes  are  identified  as  acoustic  source  terms.  The  procedure  has  produced  the  following  results: 

•  The  assumption  underlying  Lighthill’s  acoustic  analogy,  that  =  0  outside  the  source  region, 

is  justified  when  the  compressible  nature  of  the  source  term  is  included. 

•  The  compressible  strain  field  is  orthogonal  to  the  solenoidal  strain  field:  spj<p^p  =  0. 

•  A  closure  for  the  compressible  correction,  Tfjnj  =  [pUiUj  —  Poovivj]nj->  to  the  incompressible 
approximation  of  Lighthill’s  source  term,  Poo{vivj\ij  •>  has  been  obtained. 

•  Lighthill’s  acoustic  analogy  can  be  written  in  closed  form: 

C^P,tt  -  V2P  =  Poo(VzVj),ij  +  1~1qoo  [(PsViVj) nj  -2(V]~Ps),j\ 

where  — V2ps  =  poo(viVj),ij.  The  incompressible  approximation  to  the  Lighthill  source  term 
is,  of  course,  Poo{^j)nji  the  dominant  term. 

•  Though  Tg-  is  of  0(Mt2)  with  respect  to  v^j  its  contribution  to  the  acoustic  field  is  of  0(Mt) 
with  respect  to  the  contribution  from  ViVj. 
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•  The  0(M j2)  incompressible  approximation  to  the  Lighthill  source  term,  piiiUj  «  p^Vj, 
predicts  an  acoustic  field  accurate  to  O(Mt). 

The  results  indicate  the  possibility  of  investigating  the  consequences  of  compressibility  analyti¬ 
cally.  The  results  also  indicate  a  procedure  that  extends  the  possibility  of  incompressible  DNS  to 
investigate  the  compressible  effects  on  sound  generation.  This  has  some  utility  as  incompressible 
simulation  methodologies  are  better  understood,  further  advanced,  and  cheaper,  allowing  larger 
Reynolds  number  and  longer  simulation  times. 
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Appendix:  leading  order  compressible  equations 

Additional  aspects  of  the  leading  order  compressible  problem  given  in  §4  are  summarized.  The 
leading  order  compressible  momentum,  continuity  and  thermodynamic  equations  are 


D 

yr—  U>0  i  +  PO  ii 
UtQ 


=  0 


D 

Wo 


PO  +U>0  k,k 


Pa  ~  7Pc 


D 

Dh 


Ps 


-  0. 


(66) 

(67) 

(68) 


The  equations  are  nondimensional  —  the  sound  speed  is  unity.  The  compressible  velocity  gradients 
evolve  according  to: 


D 

Dta 


Wai  ,j 


—  POtij  Vpij  ) 


from  which  the  equations  for  the  divergence  and  the  vorticity  are  readily  obtained, 


(69) 


D 

-jyj~  3 '3  ~~  POtfJ  VP'>3  WQj  >P  ' 

D 

p. ,  ^0 iip  * 

Ulq 


(70) 

(71) 


On  the  slow  time  scale,  setting  jy-  =  0  one  finds  po  =  0,  v'ok  ,k  =  Ps  anf^  vpij  wo j>p  =  0.  The 
vorticity  results  solely  from  the  initial  conditions  uioq  =  0  on  the  homogeneous  solution.  Some 
wave  equations  are  derivable 


D2  _2 

^0  “  v  Po  —  vp)j  Woj,p 


(72) 


taking  jjj-  the  leading  order  equations  can  be  combined  to  give  a  third  order  wave  equation 


D  ,D2  _2  ,  D 


(73) 


For  the  slow  portion  of  the  compressible  field  w0i,j  =  0,  and  the  following  sourceless  wave 
equation  describes  the  fast  acoustic  modes: 


D 

Dto 


Po  -  V2po]  =  0. 


(74) 


In  the  general  case  ^  wq i,j  ^  0,  and  a  Lilley  type  wave  equation,  Goldstein  (1976),  results: 


D 

Dig 


D2 


\ PO  “  V  Po]  H“  Vkrt  Po  nk —  ^kn  ^p)k  ^Qiip  ■ 


(75) 


Had  the  problem  included  a  mean  shear  a  form  of  Lilley’s  equations  with  the  familiar  factors  of  2 
in  the  refractive  term  and  source  terms  would  have  resulted. 
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